Electronic properties of disordered valence-bond stripes in cuprate superconductors 
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We calculate single-particle properties of short-range ordered stripe states, using Monte-Carlo 
simulations of collective charge-density wave (CDW) order parameters coupled to fermions on a 
2d square lattice. For superconducting bond-centered stripes with a d-wave form factor, we find a 
valence-bond "glass" which coexists with low-energy quasiparticles featuring interference phenom- 
ena, in agreement with recent scanning-tunneling-microscopy (STM) measurements on underdoped 
Bi2Sr2CaCu208+i and Ca2- a: Na I Cu02Cl2. Together with earlier work, our calculations provide a 
link between CDW signatures seen in STM and those in magnetic neutron scattering. 



I. INTRODUCTION 

Charge-density wave (CDW) phenomena have been 
detected in a number of superconducting cuprates. 
Most prominent are the uni-directional spin and charge 
modulations, termed "stripes", in La2-a;Ba x Cu04 and 
La2- x Sr x Cu04 (with Nd or Eu co-doping), being 
strongest near l/8th doping. 1-5 In other cuprates, no- 
tably Bi2Sr 2 CaCu208+,5 and Ca2- s Naa;Cu02Cl2, scan- 
ning tunneling microscopy (STM) measurements have 
found signatures of short-range charge order. 6-9 . 

Remarkably, the existence of stripe states was postu- 
lated in early theory work on the Hubbard model, 10-12 
far before experimental indications for such phases were 
found. Later on, ideas of frustrated phase separa- 
tion as driving force of stripe formation were worked 
out in detail, 13 and CDW quantum criticality was pro- 
posed as source of both non-Fermi liquid behavior and 
superconductivity. 1 For the vast number of subsequent 
theoretical activities we refer to the review articles 
Ref. 15,16. 

While the role of charge order for the overall proper- 
ties of cuprates is under debate, a plausible hypothesis is 
that tendencies toward charge ordering are common to 
underdoped cuprates. 13-18 Even compounds not display- 
ing long-range order are influenced by the proximity to 
a charge-ordered state. In particular, impurities will act 
as random-field pinning centers for the collective charge 
modes, leading to static short-range order (as observed 
in STM). 15,19 ' 20 Moreover, charge order will influence the 
magnetic excitations, believed to be the pairing glue: 21 
It was recently shown 22 that short-range-ordered stripes 
give rise to an "hour-glass" magnetic spectrum, very sim- 
ilar to that observed in neutron scattering experiments 
both on La 2 - a; Ba a: Cu04 23 and YBa2Cu 3 6 +«5. 24 

While neutron and X-ray scattering were used to de- 
tect superstructure modulations from long-range charge 
order, 1,2,5 there is relatively little information on the elec- 
tronic structure of stripe states. Both STM and photoc- 
mission indicate the presence of coherent, gapless nodal 
quasiparticles (QP) in (1,1) direction, whereas antin- 
odal QP in (1,0) direction are rather incoherent and 
likely dominated by charge ordering. 8 ' 25,26 For the com- 
pound La2-xBaa;Cu04, a d- wave-like gap was recently 



reported, 2 which may be attributed to static stripes or 
to fluctuating superconductivity. 28 

In this paper, we present a detailed study of local elec- 
tronic properties of disordered stripe states in cuprates, 
using a CDW order-parameter approach plus a mean- 
field theory for the single-particle dynamics. A central 
ingredient is the <i-wave-like form factor of the charge 
order, 29 which causes the modulations to be located pri- 
marily on Cu-O-Cu bonds instead of on Cu sites. Our 
results reproduce central features of the STM data of 
Refs. 9,26. As we employ the same model for the collec- 
tive CDW modes as in Ref. 22, used there to calculate 
spin excitations in the presence of disordered stripes, our 
results provide a link between different probes of stripe 
physics. 

The remainder of the paper is organized as follows: In 
Sec. II we describe the employed model together with the 
approximations and their physical background. Sec. Ill 
presents the main numerical results, with focus on de- 
scribing the STM data of Refs. 9,26. A discussion and 
conclusion closes the paper. 



II. PHENOMENOLOGICAL MODELLING 

Our phenomenological model consists of coupled CDW 
fluctuations and electrons, with the action S — + 
S c + So/). To account for the strong commensuration 
effects observed experimentally, all fields will be defined 
for discrete lattice coordinates. 30 



A. Lattice CDW order-parameter theory 

The CDW part captures the tendency toward stripe 
ordering and is identical to that of Ref. 22: Two complex 
fields ip x ,y(^t t) represent the amplitude of horizontal and 
vertical stripe order at wavevectors K XtV , such that the 

real field Q x (f) — Reifj x (r)e' lK: ' : ' r (similarly for Q y ) mea- 
sures the modulation of both the charge density and bond 
order (i.e., kinetic energy or pairing amplitude), for r on 
sites and bonds, respectively. Then, 5p(rj) = Q x + Q y 
is the deviation of the local hole density from its spa- 
tial average. We restrict our attention to K x = (ir/2, 0) 
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and K y — (0,7r/2), i.e., a charge modulation period of 4 
lattice spacings. 6,7 The complex phase of ip x ,y represents 
the sliding degree of freedom of the density wave. 

Fluctuations of the charge order are described by a tp 4 - 
type theory S$ for the 0(4) field ip = (ip x ,ipy). The pre- 
cise form of iSy, will determine the character of the fluc- 
tuations (amplitude vs. phase). The STM data of Ref. 6, 
with modulations present everywhere in real space, point 
toward small amplitude fluctuations; in addition, the cal- 
culated spin-fluctuation spectra of Ref. 22 were only com- 
patible with experiment under the assumption of domi- 
nant phase fluctuations. Hence, we employ 



dr^2}\\d T i> lx \ 2 + \d T ip iy \ 2 + s x \ipi ; 



+ Sy\lp r: 



+C lx \lpi x ~ 1pi+x,x\ + C 2x \lpix ~ 1pi+y, X \ 
+ci y \Ay - 4>i+x.,y\ 2 + c\ v Vl>iy - lpi+ X ,y\ 2 

+ Ul^f + U 2 lpf + V\lp ix \ 2 \ipiy\ 2 

+™ (Vi+V4 4 +Vi+<)" 
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with ipi x = ip x {fi) and tp 2 = IV'izP+IV'ij/l 2 - A combination 
of u\ < and u 2 > suppresses amplitude fluctuations 
of tp. The quartic u|V'a;| 2 |V , y| 2 term regulates the repul- 
sion or attraction between horizontal and vertical stripes; 
we shall mainly employ v > leading to stripe-like or- 
der (whereas v < results in checkerboard structures). 
The phase-sensitive w term provides commensurate pin- 
ning and selects bond-centered (instead of site-centered) 
stripes 9,22 for w>0. 



Fermions 



with Q x (f i+x / 2 ) = [Q x (n) + Q x (r i+x )}/2. The coupling 
constants K1...5 decide about the electronic struture of 
the CDW state, by implementing modulations of charge 
densities and bond kinetic and pairing energies. In the 
simplest picture, stripes correspond to modulations in 
the on-site charge densities. Those are induced by k\ 
and lead to a nearly fc-independent (s-wave) CDW form 
factor (f> 2 (k) = (ct + ^ ^c^). 29,31 However, local ordering 

can instead be dominated by physics on Cu-O-Cu bonds: 
Stripe formation is driven by the competition between 
kinetic and magnetic energies, both living on bonds. 32-34 
We have recently argued 29 that such a bond-dominated 
stripe state will have modulations in (c| CT Ci+A,(r) with lo- 
cally different signs on horizontal and vertical bonds, im- 
plying a strong d-wave component of (f> 2 (k), see Fig. 1 
of Ref. 29. Modulations on bonds are induced by K2...5, 
with the d-wave character encoded, e.g., in k 2 — — K3. 

A few remarks are in order: In the advocated model, 
Eqs. (1,2,3), correlation effects are included via be- 
ing a renormalized quasiparticle dispersion and via Q XtV 
representing collective CDW tendencies, while genuine 
Mott physics is absent. Dispersion renormalizations are 
standard in mean-field theories of correlated electrons; 
here we refrain from a self-consistent calculation of the 
dispersion and instead use plausible hopping parameters 
extracted from photoemission experiments. The separa- 
tion of degrees of freedom into quasiparticlcs and col- 
lective CDW fields, both with full spatial or momentum 
dependence, is phenomenological and cannot be rigor- 
ously justified. However, e.g., in the context of electrons 
interacting with antiferromagnetic fluctuations, this has 
proven to be a fruitful route of investigation. 35 



To calculate electronic properties in the presence of 
collective charge modes, we start from a BCS model of 
fermions on the square lattice of Cu atoms: 31 

Sc = j drJ2 [% CT ( 9 r+e £ -M)c^ + ^i c kf-ki + c,c -) 



(2) 

where summation over spin indices a is implied. The 
single-particle dispersion consists of hopping to first, sec- 
ond, and third neighbors, with t — —0.15 eV, t' = —t/4, 
t" = t/\2. The chemical potential is ^ = —0.12 cV, lead- 
ing to a hole doping of « 11%. The pairing is of d-wave 
type, Ajr = Ao(cosfc K — cosfcy) with Ao = 24 mcV. 

The coupling to the collective CDW fields Q x . y reads 



+ {K2Qx(ri +x /2)CiaCi+ x ,a + K?,Q x {r i+y / 2 )c i(J Ci + y^ 
+ K4Qx(f i+x / 2 )CilC i+ xl + K 5 Q x (f i+y / 2 )c^C i+ yi+C.C. 



+ [x <-» y] 



(3) 



Observables 



STM experiments determine the spatially resolved 
local density of states (LDOS), p(f,E), up to an r- 
dependent tunnel matrix element (which depends on the 
set-point conditions. 9 ' 26 ) To separate physical modula- 
tions from set-point effects, the LDOS ratios 



Z(f,E) 
R(f,E) 



p(f,E) 
p(r, -E) ' 

fo dujpjr, uj) 

jl E dojp(f, u) 



(4) 



have been used. In a weakly doped Mott insulator, both 
Z and R (measuring spectral particle-hole asymmetry) 
can be shown to be proportional to the hole density. 36 ' 37 
In Ref. 9, spatial modulations were observed in R(r, E). 
In the following, we shall assume that these reflect mod- 
ulations in the hole density. 37,38 
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a) rf-stripe b) s-stripe c) d-checkerb d) s-checkerb 




FIG. 1: Local densities for ordered CDW states, calculated 
from S c + S c ^j for constant ip. The figures show both site and 
bond densities, (1 — c\ a c ia ) and (c| CT Ci+A,£r + /i.e.) 40 for 8 2 Cu 
sites; black squares are centers of CuC>2 plaquettes. a) Stripes 
with ki=0, K2 = — «:3 = 0.03, — «4 = K5 = 0.007 (d-wave). b) 
Stripes with ki =0.022, «2...s=0 (s-wave). c) Checkerboard 
[d- wave as in a)], d) Checkerboard [s-wave as in b)]. 



D. Perfect CDW order 

Perfectly ordered CDW states are described by S c + 
<S C .0, with ip x ,y taken to be constant. From the diagonal- 
ized fermionic Hamiltonian all electronic properties can 
be obtained. 

Sample results for the real-space densities of differ- 
ent types of CDW are displayed in Fig. 1. Here, ip x = 
(1 + i)/ v2, ip y = for the bond-centered stripes in panels 
a) and b), while ip x = i/j y = (1 + i)/V2 for the checker- 
boards in panels c) and d). The couplings k were taken 
to induce s-wave-like [panels b) and d)] or <i-wave-like 
[panels a) and c)] modulations, and the overall n ampli- 
tude was chosen such that the resulting modulation of 
fermionic densities is about 30-40%. To facilitate com- 
parison with STM data, 9 which show a strong modula- 
tion on the bonds of the Cu02 plane, 39 we included the 
bond charge densities (c| CT Ci+A,o- + h.c.) (i.e. kinetic en- 
ergies) in Fig. 1 - those are shown in between the square- 
lattice sites. 40 Clearly, Fig. la with d-wave stripes is most 
compatible with experiment. 9 



E. Pinning and adiabatic approximation 

The treatment of S c +S^+S cl p requires additional in- 
put. Pinning is important especially in the disordered 
phase of S^: Quenched disorder (e.g. from dopant impu- 
rities) acts as a random field and renders static a short- 
range ordered stripe configuration. In such a situation, 
the electronic properties can be approximately calculated 
by diagonalizing S c + S c ^p for fixed static configurations 
of ip x .y. We generate these from classical lattice Monte 
Carlo (MC) simulations of S^, , using a standard Metropo- 
lis algorithm at a finite effective temperature (T = 1) in 
a regime where the stripe correlation length is of order 
£ » 10. 41 

The numerical procedure parallels that of the adia- 
batic approximation of Ref. 22, with the difference that 
the ingredient of pinning is crucial to obtain a static sig- 
nal in STM. (The presence or absence of pinning was 



irrelevant to the finite-frequency spin fluctuations de- 
scribed in Ref. 22.) In contrast to earlier work deal- 
ing with fermionic properties in the presence of disor- 
dered stripes, 42 our modelling implements the rf-wave 
bond character, and it properly describes short-range or- 
der via Si/, (1), i.e., stripe segments coexist with checker- 
board domain walls. 19,20 ' 22 



F. Choice of parameters and validity of 
approximation 

The parameters of S c +S ci f,, Eqs. (2,3), used in our 
simulations are taken as in the static-stripe calculation 
above, i.e., for the fermionic sector we use values for t, 
t', t", and A, which are standard in the BCS mean-field 
description of cuprates, and the couplings k are taken as 
in Fig. la. 

The CDW part of the action, S$ (1), is designed 
to capture the complicated non-universal physics of the 
strongly correlated CDW formation on the lattice scale. 
The combination of s x . v , u\ and ui decides about the 
importance of amplitude vs. phase fluctuations of the 
CDW, 22 we have used s x = s y = — 4. . . — 3, u\ = —1.15, 
u 2 = 0.1. Choosing v — 0.2 prefers stripes over checker- 
boards, but allows for some checkerboard structure be- 
tween stripe domains. 22 Finally, w = 0.05 is taken for 
a moderate commensurate lattice pinning toward bond- 
centered stripes. The precise values of the mass s and the 
gradient c were used to tune the CDW correlation length 
£, which was between 10 and 30 in our simulations. Note 
that an overall scale factor in is free and determines 
the typical amplitude of ip which we have normalized to 
unity. The <S^ parameters here are identical to those used 
in Figs, la, 2a of Ref. 22 for the description of the spin 
excitations of fluctuating stripes. Moreover, the charge 
configurations generated from the MC simulations visu- 
ally match the STM results in the sense that short and 
medium stripe segments coexist with checkerboard-like 
domain walls. This property is robust with respect to 
parameter changes of 20% and more, provided that the 
correlation length £ is kept fixed. As we employ classical 
MC simulations for Sy, with time gradients absent, the 
parameters cannot be translated directly into physical 
energies or velocities. 

Our approach assumes that a mean-field picture of 
both superconductivity and charge order is a reason- 
able starting point for the description of cuprates. The 
adiabatic approximation neglects inelastic processes and 
stripe dynamics, which can be justified if the latter is 
slow (compared to the observed fermions), as happens 
in the proximity to a CDW ordering transition. Thus, 
the approximation is invalid for energies below a typical 
stripe fluctuation frequency; for strong impurity pinning 
this scale is small or zero. 41 The quasiparticle picture in 
cuprates may break down at elevated energies; as far as 
this happens due to inelastic physics, it is not captured 
by our approach (while some elastic disorder physics is 
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c) p(-0.3 eV) 1 d)Z(0.3eV) 




FIG. 2: Local observables for pinned short-range ordered 
stripes (k parameters as in Fig. la, correlation length £ ~ 10), 
showing 20 2 Cu sites of the 64 2 sample, a) Order parameter 
field Qx + Qy b) Local density (including bonds) as in Fig. 1. 
c) LDOS at -0.3 eV. d) Z(f, E) at 0.3 eV. The arrows mark 
the rows/columns along which the LDOS is shown in Fig. 5. 

captured) . 

Further, we assume that dimerization and bond order 
are the driving forces behind stripe ordering, 22 ' 34 whereas 
magnetic long-range order is less important. For simplic- 
ity, we therefore neglect both order and fluctuations in 
the triplet channel. Note that this does not mean that 
we ignore local-moment physics entirely, but instead we 
assume that those moments form singlet valence bonds, 
which is accounted for by modulated hoppings (^2,3) in 
S c ^ (3). We note that the coupling to magnetic fluctua- 
tions will contribute to the broadening primarily of antin- 
odal quasiparticlcs, 43 but a calculation including spatial 
disorder and inelastic processes is beyond the scope of 
the present paper. 




(-2x,-2x) (2jt,-2jt) 

FIG. 3: Fourier-transformed local density (from sites and 
bonds 40 ) for short-range ordered stripes, obtained from av- 
eraging 40 MC configurations on a 64 2 lattice. Left: k\ =0, 
ft 2 = -K3=0.03, _ K4 = K5 = 0.007 (dr wave). Right: ki=0.022, 
K2...5 = (s-wave). The dark/light (red/green) arrows indi- 
cate the peaks at (7r/2, 0) and (37r/2, 0). (The strong peaks at 
the zone boundary arise from artificially setting the density 
to zero in the plaquette centers.) 

Figs. 3,4 of Ref. 9. In particular, the modulation lo- 
cally breaks the C4 rotation symmetry down to C2 and 
is primarily located on the Cu-O-Cu bonds. The latter 
fact - which originates in the d-wave form factor - can 
be nicely seen in the Fourier-transformed density, Fig. 3. 
Stripe order is manifest in peaks at (ir/2, 0), (0, 7r/2) and 
(37r/2,0), (0,37r/2), with the signal at (3tt/2,0) being 
much stronger compared to (7r/2, 0) (whereas for s-wave 
stripes the peaks are roughly equal in intensity). Again, 
this is in agreement with STM data, Fig. 6 of Ref. 9. 

We note that the present comparison between theory 
and experiment does not easily allow to deduce the am- 
plitude of the actual modulations: The only observables 
free of set-point effects are Z and R. However, a reli- 
able calculation of these has to has to cope with Mott 
physics not included in our model. 38 A rough estimate, 
however, relates the experimentally observed R contrast 
of w ±30% to a bond modulation of similar magnitude. 



III. NUMERICAL RESULTS 

We now turn to the numerical results obtained from 
Sc+S^+Scip for short-range ordered stripes. 



A. Local densities 

Fig. 2 displays the order parameter field (Q x +Qy) to- 
gether with the resulting fermionic charge density, the 
LDOS p(f,E), and Z(f,E) at a high energy of 0.3 eV, 
for one fixed ip x ,y configuration 41 for a e?-wave coupling 
in S C ip. The stripe modulation, being prominent on the 
bonds, leads to a large contrast in Fig. 2b, while the 
contrast in both the site-LDOS p and Z (Figs. 2c, d) is 
weaker. (Both p and Z show strong modulations around 
the gap energy, Fig. 5 below.) The result in Fig. 2b 
has a striking similarity to the "glassy" structures in 



B. Nodal quasiparticles and quasiparticle 
interference 

One outstanding feature of the STM results on under- 
doped Ca2- x Na 2: Cu02Cl2 is the presence of quasiparticle 
interference (QPI) features in the low-energy spectra, in 
a situation where the high-energy spectra are dominated 
by period-4 modulations. 26 

Our calculations qualitatively reproduce this physics. 
For QPI to occur, we have to add realistic disorder as 
source of QP scattering: Following Ref. 44, we use a com- 
bination of 2% extended potential scatterers (strength 40 
meV, size 1.2), 5% extended pairing scatterers (strength 
Ao, size 1.5), and 0.2% pointlike unitary scatterers 
(strength 2.5 eV). The real-space results of such a cal- 
culation are displayed in Fig. 4. Both the LDOS and the 
Z map at higher energies are clearly dominated by stripe 
segments, whereas the signal below 25 meV shows 
the typical QPI modulations (compare e.g. Fig. 3c of 
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FIG. 4: LDOS p(f, E) at negative (left) and positive (middle) 
bias, together with Z(f,E) (right), for pinned short-range 
ordered stripes with additional impurities on a 64 2 lattice, at 
energies 50, 24, 8 meV (from top to bottom), k parameters 
are as in Fig. la, for details see text. 

Ref. 26). Extracting the scattering wavevectors from the 
Fourier transform of our data (not shown) is difficult due 
to the small system size; the only unambigous peak is at 
the so-called (77 wavevector (corresponding to the diago- 
nal modulations in Z at low E, Fig. 4). 

We point out two features of our results, (i) The Z 
map is more sensitive to QPI than the LDOS, because, 
to leading order, QPI modulations at positive and nega- 
tive E are anti-phase, while stripe modulations are in- 
phase. Nevertheless, the strong period-4 modulations 
seen in the experimental low-energy LDOS 26 are likely 
due to set-point effects, (ii) Real-space localization of 
antinodal QP cannot be made responsible for the loss of 
QPI at higher E. We have calculated the inverse partic- 
ipation ratio (not shown) as an indicator of localization, 
and have observed no localization signatures on scales 
up to several £ (while these length scales are sufficient to 
observe QPI). 

More generally, the compatibility of stripes with long- 
lived nodal QP has been pointed out in the past. 29 ' 34, 45-47 
For small stripe amplitude, this already follows from the 
fact that the ordering wavevector Q does not connect 
the nodal points. 34 ' 46 In our case, nodal QP survive even 
for large stripe amplitude due to the d-wave character of 
the charge order (provided that the s-wave component 
remains small). 29 

Within our simulations, the survival of coherent nodal 
QP in the presence of disordered stripes is also seen in 
the LDOS spectra in Fig. 5, taken along two different line 
cuts indicated in Fig. 2. While strong inhomogcncitics 
occur at elevated energies, in particular near the gap en- 
ergy (note the period-4 modulation in Fig. 5b around 
—50 meV), the low-energy part of the LDOS is essen- 
tially homogeneous, 47 again in striking similarity to STM 



data. 8 ' 9 (A detailed comparison of our spectra with ex- 
periment reveals several differences, which we believe to 
be related to Mott physics not captured here.) 



IV. CONCLUSIONS 

We have determined electronic properties of short- 
range ordered stripe states, coexisting with superconduc- 
tivity. Agreement with salient features of STM experi- 
ments, in particular stripy LDOS modulations at elevated 
energies coexisting with QP interference at low energies, 
is found for valence-bond stripes with d-wave-like form 
factor, singling out a specific mean-held plus stripe dis- 
order model. 

As the same collective-mode description was used ear- 
lier to model magnetic excitations in the presence of fluc- 
tuating or disordered stripes, our calculations give a uni- 
fied account of stripe signatures seen in STM and in neu- 
tron scattering, and strongly indicate that similar physics 
underlies the modulated states observed in different un- 
derdoped cuprates. 

Very recent STM experiments indicate the quasipar- 
ticle interference disappears not only at high energies, 
but at a very specific location in momentum space, ap- 
proximately at the boundary of the antiferromagnetic 
Brillouin zone. Such a feature is not part of the present 
theory, and likely requires to take into account either an- 
tiferromagnetic fluctuations or other precursors of strong 
Mott physics. 
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FIG. 5: LDOS spectra corresponding to the data in Fig. 2, 
along the two paths indicated by arrows in Figs. 2a,b. Spectra 
have been broadened by 7 meV and shifted for clarity. 
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